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Abstract 


An  analytic  study  is  made  of  the  second  order  modifications  produced  on 

the  fast  tail  electron  distribution  function  of  a  nonuniform  plasma  subjected 

to  resonant  excitation  by  wave  sources.  The  source  models  considered  can 

represent  excitation  by  external  electromagnetic  waves  propagating  obliquely 

to  the  plasma  density  gradient,  mode-conversion  of  electrostatic  whistlers, 

beat  of  two  transparent  electromagnetic  waves,  and  direct  conversion  from 

ripples  in  the  density  profile.  The  calculation  treats  the  Landau  damping 

provided  by  fast  tail  electrons  self-consistently  and  is  applicable  to  plasmas 

i o  '*/>c  y ik-iv-r  i 

having  a  long  density  scale  length  L,  i.e.,  (k^L>1/3  »  1,  where  k^'  is  the 
Debye  wave  number  of  the  warm  background  electrons.  A  threshold  condition  is 
found  for  the  formation  of  a  positive  slope  in  the  tail  distribution  by  the 
various  excitation  mechanisms.  _ 


PACS  -  52.25  Fi,  52.35  Nx,  52.40\Vlk 
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I.  Introduction 

The  modification  of  the  zero  order  electron  velocity  distribution 
function  by  resonantly  excited  waves  is  a  central  problem  encountered  in 
various  plasma  physics  applications,  e.g.,  ionospheric  modification,  laser 
fusion,  and  RF  heating.  In  some  cases  it  is  possible,  because  of  certain 
parameter  orderings,  to  decouple  the  calculation  of  the  modified  distribution 
function  from  the  determination  of  the  spatial  shape  of  the  driven  electric 
field  causing  the  modification.  Such  an  approach  is  physically  meaningful  if 
the  density  of  resonant  particles  is  small  enough  that  their  contribution  to 
the  total  energy  and  momentum  balance,  due  to  Landau  damping,  can  be  consid¬ 
ered  to  be  a  small  perturbation,  i.e.,  the  resonant  particles  behave  as  test 
particles.  An  example  of  this  type  of  non-self-consistent  calculation  has 
been  recently  reported1  by  the  authors  in  a  study  of  resonant  absorption  of 
radio  waves  in  an  ionosphere  in  which  a  small  photoelectron  population  is 
present.  In  that  earlier  study  the  shape  of  the  modulus  of  the  resonant 
electric  field  (and  hence  the  wave  spectrum)  is  predominantly  determined  by 
collisions  of  the  background  (slow)  electrons  because,  for  typical  ionospheric 
conditions,  the  effective  phase  velocity  at  the  resonance  layer  (where  the 
external  frequency  m  matches  the  local  electron  plasma  frequency  'Op(z))  is 
quite  large  compared  to  the  thermal  velocity  of  the  background  electrons. 
Hence,  the  modifications  in  the  distribution  function  can  be  calculated,  for 
large  velocities,  by  using  the  collisionally  limited  spectrum,  as  described  in 
detail  in  Ref.  1. 

In  the  present  study  we  report  a  self-consistent  calculation,  accurate  to 
second  order  in  the  electric  field  amplitude,  of  the  time-averaged  modifica¬ 
tion  produced  on  a  pre-existing  fast  electron  tail  distribution  by  resonantly 
excited  electric  fields  in  a  nonuniform  plasma.  Two  different  methods  of 


excitation  covering  a  wide  range  of  physically  realizable  configurations 
are  considered:  one  is  a  capacitor-like  pump  field  (wave  number  k0  =  Q),  and 
the  other  is  a  unidirectional  beat-type  exciter  having  finite  wave  number 
(kQ  *0).  The  k0  **  0  source  models  the  external  excitation  by  electromagnetic 
waves  propagating  obliquely  to  a  density  gradient,  as  well  as  the  internal 
excitation  by  obliquely  propagating  electrostatic  whistler  waves2,  as  may 
occur  naturally  in  the  topside  polar  ionosphere,  or  may  result  from  a  signal 
launched  from  a  spacecraft.  The  k0  *  o  source  models  the  external  pumping 
that  can  be  produced  by  the  beat  of  two  transparent  electromagnetic  waves 
having  frequencies  <i»i,  a>2  such  that  oil,  u)2  »  <*»p,  but  o»2  -  'a>i  ~  aip,  as  is 
presently  being  considered  in  laser  accelerator  schemes3.  This  type  of  source 
can  also  model  the  beat  pumping  resulting  naturally  when  an  external  electro¬ 
magnetic  wave  interacts  with  a  pre-existing  low  frequency  ripple  in  the 
density  profile,  i.e.,  the  process  known  as  direct  conversion14,5.  The  finite 
k-o  exciter  can  also  be  related  to  interactions  stimulated  by  modulated  beams, 
to  truncated  velocity  distribution  functions,  and  to  single  particle  Cerenkov 
emission . 

A  consequence  of  the  results  reported  here  is  that  the  study  of  Kef.  1 
can  be  generalized  to  conditions  in  which  a  sizeable  fast  electron  population 
(e.g.,  large  photoelectron  flux)  is  present.  In  addition,  it  is  hoped  that 
the  present  study  may  be  useful  for  undertaking  more  general  studies  (both 
theoretical  and  experimental)  of  wave-particle  interactions  in  nonuniform 
plasmas  since  a  wide  class  of  excitation  mechanisms  is  treated  self- 
consistently . 

It  is  worth  emphasizing  that  an  extensive  literature6-*1  exists  on  the 
topic  of  hot  electron  production  in  laser  fusion  experiments.  Numerous 
studies  have  demonstated  that  in  such  experiments  fast  electron  tail 


distributions  are  generated  by  intense  electrostatic  oscillations  driven  at 
plasma  resonance.  A  related  issue  which  is  presently  being  debated  by 
experts  in  this  field  is  whether  or  not  some  spectral  features  observed  in  the 
backscattered  laser  light  are  attributable 1 ^  to  distortions  in  the  fast  tail 
distribution.  Although  the  present  study  does  not  consider  the  extreme  short 
density  scale  length  environment  appropriate  for  laser  experiments,  we  report 
here  on  the  amplitude  threshold  condition  required  (according  to  second  order 
perturbation  theory)  to  create  a  positive  slope  on  the  velocity  distribution 
function  of  a  long  density  scale  length  plasma  in  which  a  zero  order  energetic 
electron  component  is  present.  We  note  further  that  various  aspects  of  this 
problem,  akin  to  some  issues  discussed  here,  have  been  recently  examined  by 
Colunga,  et  al.13,  for  the  laser  relevant  case  of  short  scale  length. 

The  physical  model  examined  in  the  present  study  consists  of  a  one¬ 
dimensional,  warm  fluid  background  plasma  having  a  zero  order  density  profile 
no(z),  which  can  be  well  approximated  in  the  neighborhood  of  the  plasma 
resonance  by  a  linear  function  of  position,  i.e.,  din  nQ(z)/dz  =«  L-1  for 
up(z)  “  **>•  The  plasma  is  subjected  to  external  pumping  by  an  effective 
electric  field  of  the  form  E0  exp[i(kQz  -  u>t )  ]  and  three  physically  different 
situations  (kQ  «  0,  kQ  >  0,  kQ  <  0)  are  investigated.  A  spatially  uniform 
population  of  energetic  electrons  is  assumed  to  be  present  with  density  n 

t 

small  compared  to  that  of  the  background  plasma,  i.e.,  nt/nG  <<  1.  The 
separation  of  the  background  scale  length  L  from  that  of  the  tail  electrons 
(i.e.,  Lt  -►  «•)  is  physically  meaningful  because  of  the  large  disparity 
between  the  collision  lengths  of  the  background  and  tail  particles.  Also,  the 
fast  particle  collision  length  is  much  larger  than  the  scale  length  of  the 
resonant  field  structures  generated  at  <*»p(z)  =  u>. 


Such  ordering  of  scale 


.%  .1  I  ft-.  1‘it‘rti.fi 


lengths  is  naturally  achieved  in  a  wide  range  of  plasma  environments,  e.g., 


laboratory  discharge  plasmas  and  the  ionosphere. 


For  illustrative  purposes,  in  the  discussion  of  the  effects  produced  by  a 


spatially  uniform  source  (kQ  =  0)  we  emphasize  two  extreme  limits  of  the 


resonantly  excited  electric  fields.  One  is  the  collisionally  dominated 


Lorentzian  shape,  aryl  the  other  is  the  convectively  limited  driven-Airy 


waveform.  Of  course,  we  also  consider  the  general  problem  that  includes 


collisions,  convection,  and  Landau  damping  by  fast  electrons.  To  obtain  a 


physical  understanding  of  the  relevant  parameter  scalings  entering  in  the 


various  cases  considered  here,  we  recall  that  the  collisionally  dominated 


plasma  resonance  exhibits  a  peak  electric  field  given  by  (w/v)  t0,  a 


characteristic  spatial  width,  (v/ui)  L,  and  an  effective  resonant  velocity 


v  *  (v/<i>)uiL,  where  v  is  the  collision  frequency  of  the  background  plasma 
L 


electrons.  In  contrast  with  this  ordering  the  convectively  limited  resonance 


(for  k0  «  o)  has  peak  electric  field  amplitude  (k^L)  E0,  characteristic  scale 


length  at  plasma  resonance  k^-1,  and  effective  particle  resonant  velocity  at 


plasma  resonance  v^  ■  w/k^,  where  kA  =  (k[)L//3)  ^ 3  L~^»  and  kp  =  oj/v  is  the 


Debye  wave  number  of  the  background  plasma  whose  thermal  velocity  is  v. 


It  then  follows  that  the  collisionally  dominated  resonance  is 


—  2/3 

operative  if  (k^L//3)  >>  w/v,  in  which  case  the  velocity  of  those  fast 


electrons  that  can  be  significantly  modified  by  the  resonance  must  be 


comparable  to  v  and  satisfies  v  /v  »  (knL//3)1/3.  Convective  limitation 
L  L 


—  2/3 

of  the  plasma  resonance  begins  to  appear  for  (k[>L//3)  ~  u>/v,  and  implies 


that  the  resonant  electron  population  that  can  be  affected  must  have  velocities 


—  — .  1/3 

v  /v  ~  (kpL//3)  «  Thus,  for  the  two  extreme  limits  of  resonant  amplitude 

A 


limitation  the  strong  wave-particle  interactions  are  anticipated  to  arise  from 


electrons  having  velocities  on  the  order  of  (kuL//3)  v,  which  can  be  large 
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for  systems  with  large  kpL,  as  is  the  case  for  ionospheric  and  space  plasmas. 


It  is  this  population  in  velocity  space  that  we  refer  to  as  the  fast  tail.  It 


is  the  goal  of  the  present  study  to  calculate  the  modifications  produced  in 


this  velocity  interval  to  second  order  in  the  amplitude  of  the  selt-consistent 


resonant  electric  fields. 


The  manuscript  is  organized  as  follows.  Section  II  describes  the 


self-consistent  calculation  of  the  resonant  electric  fields.  The  second  order 


modification  of  the  tail  distribution  and  the  threshold  for  positive  slope 


development  are  considered  in  Sec.  III.  A  summary  of  results  and  conclusions 


are  presented  in  Sec.  IV. 


II.  Self  -  Consistent  Fields 


The  geometry  of  the  one-dimensional  problem  being  considered  is  sketched 


in  Fig  1.  The  location  z  *  0  corresponds  to  the  resonant  layer  u»p  =  m  and  the 


density  scale  length  of  the  background  plasma  is  introduced  through 


2  2 

1-  <*>p  ( z ) /cu  =  z/L,  where  z  is  the  independent  variable.  In  configuration 


space  the  complex  amplitude  K(z)  of  the  self-consistent  electric  field 


driven  by  an  effective  source  field  F0  s(z)  exp(-iu)t)  (with  S(z)  representing 


the  various  pumping  schemes)  is  determined  from  Poisson's  equation,  whose 


exact  form  is  an  integral  equation  of  the  type 


/  dz'  K(z',z,u>)  E(z')  =  b0  s(z)  , 


in  which  K(z',  z,  <u)  is  a  kernel  derived  by  integrating  the  Boltzmann 


equation,  but  whose  detailed  structure  need  not  be  analysed  here.  By 


expanding  the  kernel  in  terras  of  the  ratio  of  the  background  thermal  velocity 


v  to  the  effective  phase  velocity  of  the  field  (i.e.,  the  quantity  v 


discussed  earlier),  tq .  (1)  is  approximated  in  the  lowest  order  by 
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e(z,u>)E(z)  =  E0  S(z ) 


in  which  the  local  dielectric  operator  takes  the  form 


e(z,u>)  =  ^  2  ~y  +  f  +  iT 
kjj  dz*  L 


where  T  =■  v/u>.  Equation  (3)  does  not  contain  the  Landau  damping  contribution 
of  the  background  electrons  because  for  the  case  of  interest  here,  namely 
(kDL)1/3  »  1,  it  is  negligible  since  it  scales  as  exp[-  (k[)L//3)2/3/2] . 
However,  we  are  interested  in  describing  effects  resulting  from  wave-particle 
resonances  with  a  spatially  uniform  tail  population  whose  density  is  small 
enough  that  the  location  of  the  cold  plasma  resonance  due  to  the  background 
electrons  is  not  affected.  To  include  the  Landau  damping  caused  by  the  tail 
population  we  return  to  the  general  form  of  Eq.  (2)  and  express  it  in  Fourier 
space,  namely 

id  3k 2  - 

II  diT  +  ir  "  k^r  +  i4ir  Im  xtl  E(k)  =  Eo  s(k)  » 


where 


-  oo  -ikz 

E(k)  =  /  dz  e  E(z) 


and  with  a  similar  definition  for  S(k). 

In  Eq.  (4)  the  Fourier  transform  of  the  dielectric  operator  now  contains 

the  important  imaginary  part  arising  from  the  susceptibility  of  the  tail 

electrons  x  (k).  for  a  model  Maxwellian  tail  distribution  having  the  form 
t 


fot  (V>  =  (2irvt2)1' 


exp(-v2/2v2) 


exp(-u//2k  v£)  , 
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4iam  xt(k)  =  (*gl L  )  (*) 


1/2 


0 


where  kjjt  =  kp  (v/v^)/A  is  the  tail  Debye  wave  number  and  A  =  n  /no(0).  A 
useful  feature  of  Eq.  (7)  is  that 


4wlm  Xt  00  *  gOO 

v  tt  1/2  o  o  o 

g(k)  =  A(— ^ )  kD  (  j  )  exp(-u>  /2k  )  ,  ( 

hence,  Eq.  (4)  becomes 

i  A  A  01,2  - 

Idk  E(k)  +  {i  dk  [rk  +  «(k)3  -  ^7}E(k>  =  £o  S(k)  ,  (S 


which  is  solved  by  defining  E(k)  =  exp [f(k) ]a(k)  ,  with 


F(k)- 


[fk  +  g(k) ]L 


(11 


to  obtain 

jka(k)  -  -iE0LS(k)exp[-F(k)] 


(11 


The  three  different  excitation  schemes  of  interest  here  are  contained  in  the 
choice  S(k)  =  2ir<S(k  -  kQ),  with  kG  =  0,  kQ  >  0,  k0  <  0.  He  nee ,  Eq .  (11)  can 
be  integrated  to  yield 

K(k)  =  -  2wiEoL0(k  -  k0)exp[F(k)-F(k0) ]  ,  (12 

in  which  0(k  -  k0)  refers  to  the  Heaviside  unit  step  function.  In  Eq.  (12) 
the  boundary  condition  a(-«)  =  0  has  been  chosen  because,  due  to  evanescence, 
no  waves  propagate  to  z  *  —a>. 


A.  Uniform  Pum 


Inserting  the  value  k0  =  0  in  Eq.  (12)  yields  the  electric  field  spectrum 
excited  by  a  capacitor  -  like  pump  field.  The  spectrum  is  unidirectional 
(i.e.,  k  >  0)  because  the  waves  are  excited  at  cut-otf.  In  discussing  this 
pumping  scheme  it  is  useful  to  isolate  two  extreme  limits:  in  one  the 
amplitude  is  limited  by  collisions  and  results  in  a  Lorentzian  shape  in 
configuration  space,  while  in  the  other  the  resonance  is  limited  by  convection 
and  exhibits  a  driven  -  Airy  pattern.  The  corresponding  spectra  are: 


EL(k)  =  2irEo0(k)  exp[-ij  -£—-<*  exp( -kt  2/k2 )  ]  , 


EA(k)  =  2irEo0(k)  exp[-i  j  -i  j  -  -a  exp(  -kt2/k2 )  ] 


where  k  *  (TL)-1  is  the  characteristic  wave  number  of  the  Lorentzian 
L 

resonance,  k  *  (k  L//T)2/3l-1,  a  =(w/2) 1  /2  A  k  L  (v/v  )  is  the  scaled  tail 
AD  D  t 

density,  and  k  =  w/(/2  v  )  is  the  characteristic  wave  number  of  the  tail 
t  t 

particles.  Note  that  Eq .  (13)  follows  from  Eq .  (14)  in  the  limit  (k  /k  )3  << 

L  A 

3ir/2. 

An  interesting  feature  of  Eqs.  (13)  and  (14)  is  that  | E^( k ) | 2  =  |LA(k)|2, 
hence  second  order  phase-independent  effects  caused  by  these  two  fields  are 

identical  although  their  spatial  patterns  can  be  quite  different.  In  partic¬ 
ular  the  time  -  averaged  modification  of  the  tail  distribution  function, 
accurate  to  order  | E0 |  ,  is  the  same  for  the  collisionally  and  the 
convectively  limited  resonances. 


Figure  2  exhibits  two  complementary  dependences  of  the  scaled  modulus  of 

the  spectrum  for  a  collisionally  dominated  resonance  (i.e.,  k^/kA  small). 

The  broken  curves  correspond  to  the  dependence  of  |E(k)|  on  scaled  wave  number 


(top  scale)  while  the  solid  curves  (more  useful  in  understanding  wave  -  particle 


interactions)  show  the  dependence  on  resonant  particle  velocity  (i.e.f  v  = 
w/k)  scaled  to  the  tail  velocity  vc  for  k^/kA  =  0.25  and  kt/kA  =  We  remark 
parenthetically  that  in  plotting  the  modulus  of  E(k),  kA  does  not  appear 
explicitly  but  we  scale  to  this  variable  for  convenience  in  the  latter 
analysis.  The  curves  labelled  a  =  0  are  determined  by  collisional  damping 
alone,  and  the  rest  contain  the  effect  of  self-consistent  Landau  damping.  The 
sequence  of  solid  curves  illustrates  the  role  of  increasing  tail  density  f rom _ 
a  =  0  to  50  for  fixed  kL/kA  (for  instance  if  kDL  -  104,  v/vt  -  10-1,  a  =  5u 
corresponds  to  a  tail  density  ratio  of  .04).  It  is  seen  from  Fig.  2  that  the 
principal  effect  of  self-consistency  is  to  prevent  the  acceleration  of  slower 
tail  electrons  due  to  the  quenching  of  the  large  k  components  of  the  spectrum. 
This  quenching  results  also  in  a  steepening  of  the  spectrum  as  a  function  or 

velocity,  i.e.,  3|E(k  =  m/v)|/3v  increases.  Such  an  effect  can  result  in  a 
strong  distortion  of  the  distribution  function  over  a  limited  velocity  range, 
as  is  shown  in  Sec.  III. 

The  spatial  dependence  of  the  self-consistent  fields  is  obtained  from  the 
inverse  Fourier  transform 

E(z)  -ij/'dk  E(k)e  lkz 

—00  9 

which  is  evaluated  numerically  for  the  spectra  of  Eqs.  (13)  and  (14)  to  yield 

the  results  of  Figs  3  and  4,  respectively. 

Figure  3  displays  the  spatial  pattern  associated  with  a  co 1 1 isiona 1  ly 

limited  resonance  (i.e.,  systems  with  (k  L//T)2/3  >>  u)/v).  The  a  =  0  curve 

D 

corresponds  to  E(z)/(E0L)  =  (z  +  iTL)“  ,  which  is  the  cold  plasma  resonance, 
while  the  curve  labelled  a  =  10  shows  the  effect  of  self-consistency  tor 
kL/kt  -  /3.  It  is  seen  that  Landau  damping  by  the  fast  tail  decreases  the 
peak  amplitude  of  the  collisional  resonance,  i.e.,  an  increased  effective 
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dlssipation  rate  results  whose  magnitude  can  be  approximately  estimated  from 
Eq.  (13)  by  setting  k  =  to  yield  T  >  T  {l  +  a  exp(-(TktL)2 )} .  Another 
effect  is  the  appearance  of  oscillatory  features  in  the  wings  of  the  electric 
field  modulus. 

The  spatial  pattern  of  the  convectively  limited  resonance  is  exhibited  in 
Fig.  4,  where  the  scaled  modulus  of  the  electric  field  | EA(Z ) | /kALE0  is 
plotted  versus  the  scaled  spatial  coordinate  ?  =  kAZ  =  (k.  L//3)2/3  z/L.  The 
curve  labelled  I”  =0,  a  =  0,  corresponds  to  the  undamped  mode  conversion 
result14  ir|Gi(-c)  +  i  Ai ( — C)  I  >  while  the  curve  labelled  I*’  =  0.25,  a  =  0, 
exhibits  the  effect  of  collisional  damping  and  is  given  by  n| Gi ( -?+i T ' ) |  + 
i  Ai(-£+ir')  where  Gi ,  Ai  represent  the  appropriate  Airy-type  functions15,  and 
r '  =  ^a/^L  =  (k^L//3)2/ 3( v/w)  is  the  scaled  damping.  The  bo  ttom  curve  in 
Fig.  4  contains  the  combined  effects  of  collisions  and  tail  Landau  damping  tor 
T'  »  0.25,  a  =  10,  and  (kt/kA)2  =  40.  It  is  seen  from  this  curve  that  at  this 
relatively  low  tail  particle  density  the  peak  amplitude  of  the  resonance  is 
not  affected  by  Landau  damping,  but  that  a  reduction  due  to  collisional 
dissipation  occurs.  For  the  conditions  chosen  in  Fig.  4,  tail  Landau  damping 
begins  to  overwhelm  collisional  damping  of  the  mode-converted  Langmuir  wave 
beyond  the  third  oscillation  of  the  Airy  pattern.  It  is  expected  that  at  this 
spatial  location  the  distortions  in  the  distribution  function  calculated  in 
Sec.  Ill  begin  to  appear. 

For  completeness  we  would  like  to  add  that  the  connection  between  the 
amplitude  E0  entering  into  the  kQ  =  0  model  discussed  here  and  the  power  per 


unit  area  PD  carried  by  an  electromagnetic  wave  approaching  plasma  resonance 
from  the  low  density  side  at  an  asymptotic  angle  9  is 


8nPQ  1/2 

lEo|  =  (“7T-) 


(lo) 
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i 

where  n  is  the  mode  conversion  absorption  efficiency  determined  by  Kelly  and  < 

I 

Banos1®,  and  by  Fordslund,  et  al.l/,  and  which  for  small  q  values  behaves18  as 
n  «■  a(q-2aq/3)  ,  with  a32.6  ,  and  where  q  S  (o>L/c)2/3  sin20  ,  with  c  the 
speed  of  light. 

The  corresponding  connection  between  E0  and  the  asymptotic  amplitude  tw  ■ 

of  an  obliquely  propagating  electrostatic  whistler  wave  with  frequency  u 
smaller  than  the  electron  gyrof requency  ft  ,  and  approaching  plasma  resonance 
along  Bp  from  the  high  density  side  of  a  magnetized  plasma  with  VnQ  *  B0  =  U 
is 

Eo  *  —  [cosh(tr0o)  ] 1/2  e  ^°  ,  (17) 

fH  1 

1/2 

where  0O  =  (kj_L)Y2/[  2(Y2-1 )  ]  ,  with  Y  »  fi/u>  and  k^  is  the  constant 

wave  number  perpendicular  to  the  confining  field  Bp  . 


B.  Beat  Excitation 

We  proceed  to  examine  the  resonant  excitation  of  plasma  waves  by  k0  *  o 
sources,  as  may  arise  from  various  beat  processes.  For  example,  beat  excita¬ 
tion  by  two  transparent  electromagnetic  waves  propagating  along  the  density 
gradient  in  an  unmagnetized  plasma  and  having  frequencies  u»[,  mv  ,  with  w2-u>i 
a  "^(0),  causes  a  beat-ponderomotive  force  on  the  background  electrons.  The 
value  of  the  effective  field  E0  entering  in  Eq .(13)  is  then  given  by 


E 


o 


e  <«>p(0) 
me  o»1  o»2 


Eo  1K0  2  I  exP  l  92"9  » 


(  1H) 


where  e,  m  are  the  electron  charge  and  mass,  EQj  represents  the  amplitude  ot 
the  jth  electromagnetic  wave,  and  0j  is  the  corresponding  phase.  If  the 
excitation  results  from  the  beat  of  a  long  wavelength  electromagnetic  wave  ot 

frequency  u>  »  <*^(0)  and  amplitude  ,  with  a  static  distortion  of  the 
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density  profile  having  the  form  5n  «  n  exp(ik0z 


no(0)  1 


the  equivalent  value  of  E0 

(19) 


The  magnitude  of  the  scaled  wave  spectrum  generated  by  beat  excitation  in 

the  direction  of  decreasing  density  (k0  >  q)  ±s  illustrated  in  Fig.  5.  This 

figure  shows  the  dependence  of  |E(k)|  on  both  scaled  wave  number  k/k^  (dashed 

curve)  and  on  scaled  resonant  velocity  v/v  (solid  curve)  for  k0/k  =  3, 

t  A 

T'  *  0.25,  a  ■  10,  k  /k  *  6.  In  this  case  the  fastest  phase  velocity 

t  A 

present  in  the  spectrum  is  u>/k0  (i.e.,  that  of  the  beat  source)  since  as  the 
wave  propagates  down  the  density  gradient  the  WKB  wave  number  increases  as 
k(z)  /?  .  A  consequence  of  this  dependence  is  that  the  modifications  in 
the  distribution  function  are  localized  in  velocity  space  because  no 
velocities  larger  than  <n/k0  are  present  and  as  in  the  k0  *  0  case  the  spectrum 
is  quenched  for  large  k  (slower  velocities)  by  Landau  damping.  In  comparing 
Fig.  5  with  Fig.  2,  it  should  be  noted  that  Fig.  2  corresponds  to  a 
collisionally  dominated  plasma  in  which  I"  *  4. 

The  spatial  dependence  of  the  electric  field  modulus  |E|/(k^LE0)  associ¬ 
ated  with  the  spectrum  of  Fig.  5  is  displayed  in  Fig.  6,  with  the  correspond¬ 
ing  real  part  also  shown.  In  comparing  Fig.  6  with  Fig.  4  it  is  evident  that 
beat  excitation  down  the  density  gradient  causes  a  peak  enhancement  in  the 
resonantly  excited  wave  which  is  about  a  factor  of  3  smaller  than  that 
produced  by  a  kQ  »  0  exciter  having  equal  value  of  E0.  This  reduction  has 
two  different  physical  origins:  one  is  a  larger  group  velocity  at  the 
resonance  location,  and  the  other  is  larger  Landau  damping.  The  resonance 
location  in  this  case  occurs  away  from  the  cold  plasma  resonance  layer  (i.e., 

Vo)  “  <*>)  because  the  beat  exciter  attains  resonance  in  the  neighborhood  of 


the  point  where  the  WKB  wave  number  k(z)  of  the  finite  temperature  Langmuir 
wave  matches  the  beat  kQ,  which  for  the  present  case  corresponds  to  {  ■  9, 
as  indicated  by  the  dashed  vertical  line  in  Fig.  6.  As  is  seen  from  this 
figure,  the  peak  amplitude  is  shifted  down  the  gradient  from  the  resonance 
point,  a  feature  which  is  also  present  in  the  k0  =  0  case. 

The  scaled  magnitude  of  the  spectrum  excited  by  a  beat  source  propagating 
in  the  direction  of  increasing  density  (k0  <  o)  is  shown  in  Fig.  7.  The 
parameters  used  are  the  same  as  those  of  Fig.  5,  except  for  the  sign  of  k0 f 
i.e.,  now  k0/kA  =  -3.  It  is  seen  from  Fig.  7  that  for  this  excitation  scheme 
it  is  possible  to  cause  modifications  in  the  distribution  function  of 
particles  moving  up  (v  <  0)  as  well  as  down  (v  >  0)  the  density  gradient 
because  the  resonantly  excited  wave,  initially  propagating  up  the  gradient,  is 
reflected  at  the  cut-off  point  at  z  *  0,  and  proceeds  to  propagate  down  the 
density  gradient  in  a  manner  similar  to  that  resulting  from  a  k0  =  o  exciter. 

The  spatial  dependence  of  the  scaled  electric  field  corresponding  to  the 
spectrum  of  Fig.  7  is  displayed  in  Fig.  8,  and  is  to  be  compared  with  the 
behavior  shown  in  Fig.  6. 
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III.  Modification  of  Tall  Distribution 

We  proceed  to  calculate  the  modification  produced  in  the  pre-existing, 
spatially  uniform,  fast  tail  electron  distribution  by  the  various  self- 
consistent  fields  investigated  in  Sec.  II.  The  quantity  of  interest  here  is 
the  time-averaged,  spatially  asymptotic  change  in  the  distribution  function, 
accurate  to  second  order  in  E0.  The  meaning  of  the  time  average  refers  to 
time  samples  long  compared  to  the  period  of  the  coherent  oscillation,  i.e., 

At  »  2h/(d0.  Spatially  asymptotic  implies  that  the  region  of  interest  is  far 
from  the  spatial  point  where  the  wave  particle  resonance  is  encountered  for  a 
given  velocity,  i.e.,  where  k(z)  =  <o/v  .  The  restriction  to  second  order 
perturbation  theory  implies  that  strong  nonlinear  wave-particle  interactions, 
such  as  particle  trapping,  are  not  included. 

It  is  useful  to  express  the  time-averaged  tail  distribution  function  in 
the  form 

±  + 

<  fc  (z  »  ±«,v)  >  =  fot(v)  +  <  6f ^  (z  =  ±<*>,v)  >  ,  (20) 

where  the  ~  label  refers  to  v  >  0  and  v  <  0  particles,  respectively,  and  the 

<>  represents  the  time-averaging  operation. 

The  modification  in  the  distribution  function  is  obtained  by  integrating 

the  time-averaged  Vlasov  equation  in  configuration  space  using  appropriate 

+ 

boundary  conditions  for  v  >  0  and  v  <  0,  namely  <6ft  (z  =  -<»fV)>  = 

<6ft  (z*«*,v)>  »  0.  The  relevant  expression  for  v  >  0  is 

00 

<  6f£  (z~,v)>  -  «  {/  dz'E(z')  [ft  (z',v)]*  }  +  c.c.  ,  (21) 

4mv  dv  i 

—00 

where,  f*  (z',v)  is  the  first  order  modified  distribution  function  oscillating 
as  exp(-iwt)  and  given  by 
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a  z 

(z’.v)  =*  —  —  fQc(v)  /  dz"  exp  [l  —  (z  ’-z")  ]  £(z")  .  (22) 

1  mv  av  „  v 

Inserting  Eq .  (22)  into  Eq.  (21)  and  representing  the  spatially  dependent 
electric  fields  in  terms  of  the  Fourier  transform  defined  by  Eq.  (15)  yields 
the  familiar  expression 


< &ft  ( z 5 


<li(k  •r)lIit'o«J. 


v  3v 


(23) 


in  which  it  has  been  consistently  assumed  in  all  intermediate  integrations 
that  u  has  an  infinitesimally  small  imaginary  part  that  causes  the  oscillatory 
contribution  at  z  =*  to  vanish.  A  similar  procedure  using  the  appropriate 
boundary  conditions  for  v  <  0  particles  yields  an  identical  expression,  thus 
permitting  the  generalization 


<«ft-  (z=±-,v)>  =  (y  i 


3_ 
v  3v 


I JL 

v  3v 


cot 


(2A) 


It  should  be  mentioned  that  the  expression  in  Eq .  (24)  is  often  referred  to  as 
the  quasilinear  result,  but  in  the  present  context  it  should  be  properly 
interpreted  as  the  spatially  asymptotic,  second  order,  time-averaged  modifica¬ 
tion  produced  by  a  temporally  coherent  exciter. 

A.  Uniform  Pump 

The  spectrum  excited  by  a  k0  =  (j  source  in  a  nonuniform  plasma  is  strictly 
unidirectional,  as  shown  in  Eqs.  (13)  and  (14).  Consequently,  the  modifica¬ 
tions  in  the  tail  distribution  function  are  restricted  to  those  particles  that 
propagate  in  the  direction  of  decreasing  density,  which  implies  that  <Sf->  =»  0 
or 

<ft  (z— ,v)>  -  f ot ( v )  ,  (25) 


i 
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and 


£  2 

<f£(z««°,v)>  -  fot(v)  [l  +  )  h+(v )  ]  , 


(2b) 


where 


vt2vL 


h+(v)  =  {  1 — 2[ — 5 —  +  a  exp(-v2/2v2 )  ] }  * 


exp{-2f^-  +  a  exp(-v2/ 2v^  )  ] }  , 


2/  2 ' 


(27) 


and  in  which  the  quantity  Es  =  mv^2/(effL)  represents  the  natural  scaling  for 


the  strength  of  the  pump  field  EQ  that  can  cause  significant  tail  modifica¬ 
tions.  In  the  manipulations  leading  to  Eqs.  (2b)  and  (27)  the  explicit  form 


of  t  (v)  given  by  Eq.  (6)  has  been  used  and  we  note  that  v  =  TuiL  =»  u>/k  . 
ot  L  L 


It  should  be  emphasized  that  Eqs.  (26)  and  (27)  apply  to  both  extreme 
limits  of  resonant  exicitation  by  a  kQ  =  0  pump,  i.e.,  to  the  Lorentzian  and 
to  the  driven  Airy  pattern.  In  spite  of  the  Important  role  played  by  the 
characteristic  wave  number  kA  tn  determining  the  spatial  pattern  seen  in  the 
curves  of  Fig.  4,  it  plays  no  role  in  Eqs.  (2b)  and  (27),  i.e.,  the  modifica¬ 
tions  in  the  distribution  function  are  independent  of  kt/kA  .  The  relevant 
scaled  parameters  that  determine  second  order  modifications  of  the  tail 


distribution  by  a  kQ  =  0  exciter  are:  E0/Es,  v  /v  ,  and  a. 

L  t 


In  the  large  velocity  limit,  i.e.,  v  »  v  and  v  >>  /2  v  , 

L  t 


-a.  £  2 

<ft+  (z— ,v)>  ♦  fOC(v)  (l  +  [-g2-]  ] 


(28) 


which  implies  an  enhancement  of  the  density  of  extremely  energetic  electrons 
with  no  change  in  the  shape  of  their  velocity  distribution.  This  behavior 
occurs  because  the  derivative  of  the  modulus  of  the  spectrum  witli  respect  to 
velocity  goes  to  zero  at  large  velocities  and  the  velocity  derivative  ot  a 


Maxwellian  tail  distribution  goes  as  3f  / 3v  +  -vf  /v  2  so  that  the  spatially 

ot  ot  t 


asymptotic,  time-averaged  change  has  the  form 


with  D  =  (e/2m)2  |  E(k*u)/v)  |  2  nearly  constant,  hence 


<6it  (z  — ,v)>  *  —  fQt(v) 


which  indicates  an  energetic  tail  density  enhancement  of  n  D/v  14 ,  independent 

t  t 

of  v  and  a  . 

la 

Figure  9  illustrates  various  properties  of  the  second  order  modifications 
produced  by  a  kQ  ■  0  exciter;  it  displays  logarithmic  plots  of  the  spatially 
asymptotic,  time-averaged  distribution  function  versus  scaled  velocity 
v//2  v  for  different  conditions.  The  curve  labelled  EQ  «  o  is  the  zero 
order  tail  distribution  and  the  solid  curve  marked  a  *  U  provides  an  example 
of  a  non-self-consistent  calculation  (spectrum  determined  entirely  by  colli¬ 
sions)  for  v  //2  v  “  3,  (E0/Es)2  a  35.  The  dashed  curve  is  a  self-consistent 
L  t 

result  for  a  scaled  tail  density  a  =  30.  It  is  seen  from  Fig.  9  that  tor  large 
velocities  the  self-consistent  and  non-self-consistent  results  are  identical 
and  cause  a  direct  enhancement  in  the  energetic  tail  density,  as  explained 
by  Eq .  (30).  The  modification  produced  by  the  non-self-consistent  spectrum  is 
gentler  in  slope  and  extends  to  lower  velocities  than  that  caused  by  the 
self-consistent  spectrum,  as  anticipated  from  the  behavior  displayed  in  Fig. 

-  .2 

2.  The  increased  value  of  3|E(k=u>/v)|  /  3v  due  to  self-consistent  depletion 
of  the  large  k  components  of  the  spectrum  by  Landau  damping  gives  rise  to  a 

strong  modification  of  the  velocity  distribution  over  a  narrow  velocity 

interval.  In  particular,  as  the  value  of  the  parameter  a  is  increased,  for 


fixed  Eo/Es,  it  is  possible  for  the  tail  distribution  to  develop  a  positive 
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slope  over  a  restricted  velocity  range,  as  is  illustrated  in  Fig.  10  for 


(Eq/Es)2  *  60*  v  //2  v  -3.  It  is  seen  from  Fig.  10  that  the  non-self- 


consistent  calculation  (solid  curve  labelled  a  *  0)  predicts  a  monotonic 


distribution  function,  while  the  self-consistent  analysis  yields  a  burap-on- 


tail  distortion  for  a  *  50.  For  reference,  the  unperturbed  distribution  is 


indicated  by  the  dashed  curve  (E0  =  0). 


Since  the  requirements  for  slope  reversal  with  a  kQ  =  0  exciter  have  been 


discussed  in  our  earlier1  non-self-consistent  study  in  which  the  modulus  of 


the  spectrum  is  determined  entirely  by  collisions,  we  want  to  emphasize  that, 


for  a  self-consistent  field  in  a  collisionless  plasma  slope  reversal  can  also 


result  within  second  order  theory,  as  is  illustrated  in  Fig.  11  for  T  =  0, 


a  -  10,  and  a  relatively  small  pump  amplitude  E0/Es  =  /5. 


The  threshold  condition  on  the  pump  amplitude  E0  required  to  develop 


slope  reversal  is  obtained  in  the  general  case  by  setting 


?  <f*  (z”**»v)>  =  0 


and  using  Eq .  (26)  to  solve  for  (E0/Es)2 


(?)  ■JT[fh+(v)-^hnv>r 

fc-s  V,-2  3v  V  2 


which  using  Eq.  (27)  takes  the  form 


(?)  *  b  exp((2vL/vt)  +  2a  exp( -v  2/ 2vt  2) 


where 


At  it  **•  < 


..i  'itk  •«!  i 


1  a)  '»l. -at  1  af  - 1|  -  »tl 


b  =  (6vt4v  /v5)  +  2a  exp(-v2/2vt2)  - 

Li 

2 

[l  -  (2vt2v  /v2)  -  2a  exp(-v2/2vt2)  J  .  (34) 

Li 

Since  the  left-hand  side  as  well  as  the  right-hand  exponential  factor  in  Eq. 

(33)  are  positive  definite,  it  follows  that  slope  reversal  within  second  order 

perturbation  theory  can  occur  only  for  that  restricted  velocity  interval  in 

the  distribution  function  over  which  the  quantity  b  given  by  Eq .  (34)  is 

positive.  In  general  this  range  in  v  depends  on  the  parameters  v  /v  and  a 

L  t 

through  a  transcendental  equation  that  must  be  solved  numerically.  However, 
in  the  non-self-consistent  limit  (a  *  0),  the  range  satisfies  the  inequality 

ilk  .  <  (V_)2  <  irt  ,  ,  (35) 

v  **  v  V  V  v  V 


while  in  the  collisionless  limit  (v  =0)  b  is  positive  in  the  range, 

L 

-  2  _ 

In  [ 4a/ ( 3  +  /5)]  <  — )  <  In  [4a/(3-  /5)  ] 

/2  v. 


Because  the  inequality  of  Eq.  (36)  exhibits  a  logarithmic  dependence  on  a  the 
range  over  which  slope  reversal  can  develop  in  the  collisionless  limit  is 
narrower  than  in  the  non-self-consistent  (a  =  0)  limit.  Thus  a  threshold 
field  exists  for  producing  slope  reversal  but  only  over  a  limited  range  in 
velocity  space. 


B.  Beat  Excitation 

For  a  beat  source  having  k0  >  0  the  resulting  electric  field  spectrum 
given  by  Eq.  (12)  is  unidirectional,  as  is  the  case  for  a  k0  =  0  exciter,  but 
the  fastest  Fourier  component  present  has  phase  velocity  v0  =  u)/k0. 
Consequently,  the  modifications  in  the  fast  electron  distribution  function 
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caused  by  this  excitation  process  are  restricted  to  electrons  moving  in  the 
direction  of  decreasing  density  within  the  interval  0  <  v  <  vQ  ,  i.e., 

<fc  (z*-*,v ) >  =  fot(v)  ,  (39) 

and 


<f^  (z=*®,v)> 


=  fot(v)  [  1  + 


h+(v)0(vo 


v>; 


(40) 


where  h+(v)  is 
from  Eq.  (12) 
E0  in  the  k0  = 

E+  . 


given  by  Eq.(28),  8  is  the  Heaviside 

the  tactor  exp(-F(k0)]  enters  and  the 

0  case  is  replaced  by 
v^ 

Eo  exp[~ ~  +  «  exp(-v02/2vt2 ) ) 
vo 


step  function,  and  now 
role  of  the  pump  strength 


(41) 


A  significant  difference  from  the  effect  caused  by  a  k0  =  0  pump  is  that 
the  enhancement  in  the  density  of  extremely  energetic  electrons,  given  by 
Eq .  (28),  is  not  present.  In  addition,  the  formation  of  a  positive  slope  in 
the  distribution  function  requires  the  phase  velocity  of  the  beat  source  to  be 
larger  than  the  smallest  velocity  for  which  the  factor  b  of  Eq .  (34)  vanishes, 
e.g.,  in  the  col lisionless  limit  it  requires  that 

C— 2 — )2  >  In  [  4a/ (  3+/ 5)  ]  .  (42) 

/2  vt 

For  beat  sources  satisfying  this  condition  the  threshold  electric  field  for 

bump  formation  is  given  by  the  expression  in  Eq.  (33),  but  with  E0  replaced  by 

E+.  The  characteristic  features  of  modifications  produced  by  a  kQ  >  0  beat 

source  are  summarized  in  Fig.  12  for  v0//2  v  =  3.3  ,  v  // 2  v  =  4  ,  a  =  10, 

t  L  t 

and  ( Eq/Es ) 2  «  30. 

A  beat  source  having  kQ  <  u  gives  rise  to  a  bi-directional  electric 
field  spectrum,  but  its  effect  on  the  tail  distribution  function  is  asymmetric 


with  respect  to  particles  moving  up  or  down  the  density  profile.  For 
particles  with  v  >  0  the  modification  is  similar  to  that  caused  by  a  kQ  =  0 
exciter,  i.e.. 


<*t  (z-«,v)>  =  f  ot  ( V )  [  1  +  (1^)  h+(  V )  ]  ,  (4l) 

cs 

with  the  pump  strength  E0  now  replaced  by 

V 

L  2 

E-  -  E0  exp[-  |Vq|  +  a  exp{-v02/2vt )  ]  .  (42) 

For  v  <  0  particles  the  modification  is  restricted  to  velocities  having 
magnitude  larger  than  |vQ|  t  because  the  resonantly  excited  wave  initially 
having  phase  velocity  -u/jk0j  =  vQ  approaches  cut-off  and  the  phase  velocity 
increases,  i.e,  in  the  WKB  sense  o>/k(z)  ♦  -*  as  z  +  0+.  The  modified  distrib¬ 
ution  function  in  the  range  v  <  G  is 

<f~  (z— ,v)>  =-  f0t(v)  [l  +  (|^)  h+(v)8(|v|  -  |vQ  | )  ]  .  (43) 

In  this  excitation  scheme  there  is  an  identical  energetic  electron  density 

enhancement  n^ [  1  +  (E_/es)2]  produced  up  and  down  the  density  profile.  For 

v  >  0  particles  the  condition  for  bump  formation  is  similar  to  that  of  the 

k0  *  G  pump  except  for  the  replacement  of  E0  by  E-  in  Eq.(33).  An  example  of 

a  distribution  modified  by  a  kQ  <  Q  exciter  is  given  in  Fig.  13  for  vQ// 2  vj;  = 

-3,  v  //2  v  -  10~3  ,  a  *  1G,  (E0/Es)2  =  /6.  Note  that  in  Fig.  13  the 
L  t  ± 

distribution  functions  <f  >  are  presented  for  illustrative  purposes  as  a 
single  function  of  velocity,  but  it  should  be  understood  that  the  domain  v  <  0 


is  found  spatially  at  z  +  -< 


while  v  >  0  corresponds  to  z  >  a>. 
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IV.  Conclusions 

An  analytic  study  has  been  made  of  the  self-consistent  modifications 
produced  on  a  pre-existing  fast  electron  distribution  function  in  a  nonuniform 
plasma  which  is  subjected  to  resonant  excitation  by  external  wave  sources. 
Since  the  wave-particle  interaction  analyzed  here  is  restricted  to  a  fast  tail 
distribution  which  is  initially  spatially  uniform  while  the  nonuniform  warm 
background  plasma  undergoes  collisional  heating,  the  results  are  applicable  to 
experiments  in  which  the  plasma  has  a  long  density  scale  length,  i.e., 

(kj)L)1/^  >>  1  .  Such  a  parameter  ordering  is  realizable  in  the  laboratory, 
and  is  naturally  achieved  in  ionospheric  and  space  plasmas. 

Although  the  time-averaged  modification  in  the  distribution  function 
has  been  determined  self-consistently ,  the  calculation  has  been  performed 
only  to  second  order  accuracy  in  the  strength  of  the  electric  fields.  This 
implies  that  strongly  nonlinear  effects  such  as  particle  trapping  and 
nonlinear  shifts  in  the  wave-particle  resonance  condition  v  =  ui/k(z)  are  not 
included. 

One  of  the  interesting  results  of  the  present  analysis  is  that  a  positive 

slope  can  be  formed  in  the  distribution  function  by  self-consistent  second 

order  nonlinearities  in  a  nonuniform  plasma.  The  cause  of  this  effect  can  be 

traced  to  the  competition  between  the  negative  slope  of  the  zero  order 

distribution  function  f  (v)  and  the  positive  slope  of  the  self-consistent 

ot 

spectrum,  i.e.,  3  |  E(k=u)/v )  |  2/ 3v  >  0  .  Physically,  these  opposing  tendencies 

result  in  the  existence  of  a  restricted  velocity  interval  within  which  more 
particles  are  accelerated  out  of  the  Interval  toward  larger  velocities  than 
enter  into  it  from  lower  velocities,  provided  the  pump  field  amplitude  E0 
is  above  a  threshold  value.  It  is  found  here  that  self-consistency  increases 
the  value  of  3 1  E(k=u>/v )  |  2/ 3v  and  hence  lowers  the  threshold  tor  bump  forma- 
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Cion  from  the  value  predicted  by  a  non-self-consistent  analysis.  This  feature 
is  worth  emphasizing  because  on  first  thought  one  might  expect  the  opposite 
behavior,  i.e.,  self-consistency  reduces  the  amplitude  of  the  spectrum  and 
hence  the  modifications  should  be  weaker.  That  is  not  the  case  because  the 
quantity  of  relevance  is  the  slope  of  the  spectrum,  which  is  found  to 
increase  as  the  scaled  tail  density  parameter  cr  =  (nc/n0 )( ir/2) 1  /2( v/v^ )  kyL 
increases . 

A  significant  experimental  consequence  of  slope  reversal  in  the  fast  tail 
distribution  function  is  the  spontaneous  growth  of  waves  (or  sidebands)  having 
frequencies  different  from  that  of  the  monochromatic  external  signal  driving 
the  linear  plasma  resonance.  It  is  important  to  assess  the  role  of  such  an 
effect  when  interpreting  experimental  observations  of  extraneous  frequencies 
in  resonant  absorption  experiments  in  plasmas  since  its  presence  can  be 
confused  with  parametric  instabilities.  While  in  numerous  cases  parametric 
instabilities  are  indeed  present,  it  is  possible  that  some  of  the  observed 
spectral  features  actually  arise  from  wave  induced  distortions  in  the  distrib¬ 
ution  function.  As  suggested  by  the  present  study  this  is  a  topic  that 
deserves  closer  scrutiny  in  the  context  of  specific  experiments.  In  fact, 
experiments  * ^ ^ ®  performed  in  uniform  plasmas  excited  with  grid  antennas  have 
been  found  to  exhibit  an  analogous  effect. 

In  the  study  of  the  self-consistent  resonant  electric  fields  two  possible 
excitation  models  have  been  considered  which  describe  a  broad  class  of 
experimentally  interesting  scenarios.  For  instance,  the  k0  =  0  exciter  can 
model  the  irradiation  of  an  unmagnetized  plasma  with  an  electromagnetic  wave 
propagating  obliquely  to  the  density  gradient.  This  same  model  also  describes 
the  fields  produced  during  mode  conversion  of  an  oblique  electrostatic 


whistler  into  a  Langmuir  wave  in  a  magnetized  plasma.  The  other  excitation 
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mechanisms  considered  (kQ  >  0,  k0  <  0)  model  sources  which  experimentally  can 
arise  from  the  beat  of  two  transparent  electromagnetic  waves,  as  is  presently 
being  considered  in  the  context  of  laser  accelerators,  or  may  result  from  the 
beat  of  an  external  electromagnetic  wave  with  low  frequency  distortions  of  the 
density  profile.  For  completeness,  in  the  present  work  we  have  also  indicated 
how  the  effective  pump  strength  E0  related  to  the  various  physical 
parameters  for  these  different  resonant  excitation  schemes. 

Beat  sources  that  propagate  in  the  direction  of  decreasing  density  with 
velocity  vQ  >  0  are  found  to  cause  distortions  in  the  tail  distribution 
within  the  interval  vq  >  v  >  0  because  the  WKB  phase  velocity  of  the 
resonantly  excited  wave  decreases  as  it  propagates  down  the  density  gradient. 
Although  it  is  possible  to  obtain  second  order  slope  reversal  in  this  scheme, 
the  possibility  is  restricted  by  the  value  of  vQ.  The  effect  of  a  beat 
source  propagating  in  the  direction  of  increasing  density  v0  <  0  is  highly 
nonlocal.  For  z  +  -«•  it  causes  a  modification  in  the  distribution  function 
for  fast  particles  satisfying  v  <  v0  <  o  ,  while  simultaneously  modifying  the 
distribution  function  over  the  entire  velocity  range  at  z  >  °°. 

In  summary,  various  features  of  self-consistent  wave-particle  inter¬ 
actions  in  nonuniform  plasmas  have  been  ellucidated  by  this  analytic  study. 

It  is  hoped  that  some  of  the  results  will  help  in  the  interpretation  of 
laboratory  and  ionospheric  experiments  and  stimulate  more  advanced  analysis 
of  kinetic  processes  in  nonuniform  media. 
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vW 


Fig.  1 


Schematic  of  geometry  considered.  The  background  density  profile 


Fig.  2 


Fig.  3 


Fig.  4 


nQ(z)  is  approximately  linear  near  the  plasma  resonance  located 
at  z«  0,  while  the  fast  electron  density  profile  nt  is  uniform. 
The  characteristic  velocity  of  tail  electrons  scales  as 

vt  “■  (kpL)1/3  v  and  the  various  excitation  scenarios  are 
modelled  by  S(z). 

Modulus  of  spectrum  excited  by  kQ  =  o  source.  Dashed  curves  show 

wave  number  dependence  (top  scale)  and  solid  curves  resonant 

velocity  dependence  (bottom  scale).  a  =0  corresponds  to  colli- 

sional  damping  alone,  and  finite  a  indicates  effect  of  self- 

consistency  for  k  /k  »  0.25,  k  /k  «  6  . 

Li  A  t  A 

Spatial  dependence  of  scaled  modulus  of  collisionally  limited 
resonant  electric  field  excited  by  k0  =  o  source.  Solid  curve 
shows  dependence  in  the  absence  of  fast  electrons  and  dashed 
curve  exhibits  the  self-consistent  effect  of  tail  Landau  damping 
for  a  scaled  tail  density  a  =  10  and  k t /Vcl  =  /3  .  kL  =  (co/v)L_1 
and  a  is  defined  after  Eq.  (14). 

Spatial  dependence  of  scaled  modulus  of  electric  field  for  a 
convectively  limited  resonance.  Top  curve  is  the  undamped  mode 
conversion  result  x|Gi(-5)  +  i  Ai(-c)|»  and  the  curve  I"  =  0.25, 
a  =*  10  shows  the  effect  of  collisions.  The  bottom  curve  contains 
collisional  and  fast  electron  tail  Landau  damping  for  I"  =  0.25, 
a  =  10,  and  (kt/k^)2  =  40  •  k^  5  (k[)L//3)  2/3  L-1,  and  a  is  the 
scaled  tail  density  defined  after  Eq .  (14). 
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Fig.  5 


Fig.  6 


Fig.  7 


Fig.  8 


Fig.  9 


Modulus  of  spectrum  excited  by  a  beat  source  propagating  in  the 


direction  of  decreasing  density  with  k  /k  =  3  .  Dashed  curve 

o  A 


(top  scale)  shows  dependence  on  scaled  wave  number  and  solid 
curve  (bottom  scale)  the  resonant  velocity  dependence  k  =  u>/v 


for  I”  »  0.25,  a  =*  10,  k  /k  =  6  . 

t  A 


Spatial  dependence  of  scaled  modulus  and  real  part  of  electric 
field  excited  by  a  beat  source  propagating  in  the  direction  of 
decreasing  density  with  k0/kA  =  3  and  corresponding  to  the 
spectrum  shown  in  Fig.  5.  The  dashed  vertical  line  indicates  the 
WKB  wave  resonance  point,  i.e.,  where  k(z)  =  kQ  . 


Scaled  magnitude  of  spectrum  excited  by  a  beat  source  propagating 
in  the  direction  of  increasing  density,  k0/kA  *  -3  ,  T'  -  0.25  , 


k  /k  =  6,  a  =  10  .  Dashed  curve  (top  scale)  is  the  wave  number 
T  A 


dependence  and  solid  curve  (bottom  scale)  the  resonant  velocity 
dependence  k  =  oi/v  . 


Spatial  dependence  of  scaled  modulus  and  real  part  (dashed  curve) 
of  electric  field  excited  by  a  beat  source  propagating  in  the 
direction  of  increasing  density  with  k0/kA  =  -3  and  corresponding 
to  the  spectrum  shown  in  Fig.  7.  The  dashed  vertical  line  indi¬ 
cates  the  WKB  wave  resonance  point,  i.e.,  where  k(z)  =  k0  . 


Distortion  in  tail  velocity  distribution  function  caused  by  a 
kQ  «*  0  exciter.  The  curve  labelled  E0  =  0  is  the  zero  order  dis¬ 
tribution,  while  the  a  =  0  gives  the  non-self-consistent  result. 
The  dashed  curve  is  a  self-consistent  calculation  for  a  =  30, 


v  //2  v  =  3  ,  (Eo/Es)2  _  35. 

L  t 


■  i 

V, 


y 

*  . • .  *V» 3 


Fig.  10  Example  of  localized  distortion  of  distribution  function  caused 

by  self-consistency.  Dashed  curve  is  zero  order  distribution, 

curve  labelled  a  =  0  is  non-self-consistent  result  for 

v  //2  v  =3,  (E0/Es)2  =  60,  and  curve  labelled  a  =  50  is 
L  t 

self-consistent  calculation. 

Fig.  11  Example  of  collisionless  (T  =  0)  slope  reversal  caused  by 

relatively  small  pump  field  Eq/es  =  /5.  Dashed  curve  is  the 
zero  order  tail  distribution  and  solid  curve  the  modified 
distribution  for  a  scaled  tail  density  a  =  10  . 

Fig.  12  Modification  of  the  tail  distribution  function  produced  by  a  beat 

source  propagating  in  the  direction  of  decreasing  density  with 

v0//2  v  *  3.5  for  v  // 2  v  =  4  ,  a  =  10,  (E0/Es)2  =  30.  The 
t  L.  t 

dashed  curve  is  the  zero  order  distribution. 

Fig.  13  Modifications  of  the  tail  distribution  function  produced  by  a 

beat  source  propagating  in  the  direction  of  increasing  density 

with  vq//2  v  =  -3  for  a  =  10,  v  //2  v  =  10”^,  (E0/Es)  =  /b. 
t  L<  t 

Note  that  the  modifications  for  v  <  0  appear  as  z  ♦  -°»,  while 
those  for  v  >  0  correspond  to  z  >  ".  The  dashed  curve  is  the 


zero  order  tail  distribution 


I  t  »  *  Ij  ♦  t.t  ■ 


.«  «.»  i.i  .1  ' 


T3^Z 
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